High-Order Harmonic Generation and Molecular Orbital Tomography: 
Characteristics of Molecular Recollision Electronic Wave Packets 
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We investigate the orientation dependence of molecular high-order harmonic generation (HHG) 
both numerically and analytically. We show that the molecular recollision electronic wave packets 
(REWPs) in the HHG are closely related to the ionization potential as well as the particular orbital 
from which it ionized. As a result, the spectral amplitude of the molecular REWP can be significantly 
different from its reference atom (i.e., with the same ionization potential as the molecule under 
study) in some energy regions due to the interference between the atomic cores of the molecules. 
This finding is important for molecular orbital tomography using HHG [Nature 432, 867(2004)]. 
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The orientation dependence of molecular high-order 
harmonic generation (HHG) in strong laser fields con- 
stantly attracts attention in both experiment and 
theory [H H 0, i, 0, @, 0, S] . A new surge was brought on 
by the recent emergence of its important application in 
molecular orbital tomography. Itatani et al showed that 
through calibrating the recollision electronic wave packet 
(REWP) by a reference atom, the orientation dependence 
can be used to reconstruct the shape of the highest oc- 
cupied molecular orbital (HOMO) [9|] . The molecular or- 
bital tomography is based on the assumptions that, i) 
a reference atom exists that is mainly characterized by 
an ionization potential close to that of the molecule and, 
ii) the spectral amplitude of the REWP for the reference 
atom is always identical to that of the molecule, indepen- 
dent of the alignment of the molecule [H, [l(| EL fl2l Il3|. 
Investigation of the foundations of molecular orbital to- 
mography has led to valuable insights into the mechanism 
of atomic and molecular HHG. For example, Levesque 
et al.pdj conjectures a universality in the behavior of 
atomic HHG spectra scaled by the recombination tran- 
sition atomic dipole moment. This conjecture and many 
other assumptions connected to the idea of molecular or- 
bital tomography call for a thorough theoretical exami- 
nation. 

In the present paper, we directly examine the second 
assumption by numerically calculating the spectral am- 
plitudes for ID and 2D H^ molecular ions with differ- 
ent internuclear distances and laser intensities. Through 
comparison with reference atoms that share the same 
ionization potential, we find that the REWP of Hj is 
strongly dependent on its orientation in certain energy 
regimes where the interference between the two atomic 
cores of the molecules dominates 15| . Analytically, we 
apply the Lewenstein theory to N2 and CO2 molecules 
that have different HOMOs and internuclear distances. 
Our theoretical calculations validate the assumption for 
N2 in the plateau region of the harmonic spectra that 



are of most crucial for tomographic reconstruction of the 
orbital. But, our calculations invalidate the assumption 
for CO2, showing a strong orientation dependence of the 
spectra in the plateau region due to the interference. This 
finding suggests that for molecules such as CO2, to accu- 
rately reconstruct its orbital using the HHG, the inter- 
ference effects need to be considered. 

The Hamiltonian of H^ or hydrogen-like atoms studied 
here is H(t) = p 2 /2 + F(r) - r • E(t), with soft-Coulomb 
potential V(x)= , — ~ z H . . ~ z and 
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and 2D case, respectively, where Z is the effective charge, 
R is the internuclear separation for H^", and R = a.u. 
corresponds to the hydrogen-like atom, (for Z = 1 and 
R = 2 a.u., the ground state energy for H^ reproduced 
here is Eo = 1.11 a.u. in both cases). E(i) is the exter- 
nal electric field. Here, we suppose that the molecular 
axis is coincident with the x-axis and the external field 
is linearly polarized with an orientation angle of 9 to the 
molecular axis. In the following calculations, the atom 
units of h = e = m e = 1 are adopted. Our calcula- 
tion will be performed for 780 nm trapezoidally shaped 
laser pulses with a total duration of 10 optical cycles and 
linear ramps of three optical cycles. Numerically, the 
above Schrodinger equation is solved by the operator- 
splitting method. The fast Fourier transform (FFT) is 
employed to transform the wave function between the 
position representation and the momentum representa- 
tion. The boundary is set as ±200 a.u. and the grid 
number is 4096. The ground state is obtained by propa- 
gation in imaginary time. The numerical convergence is 
checked with using a finer grid. The coherent part of the 
harmonic spectrum is obtained from the Fourier trans- 
formed dipole acceleration expectation value, and only 
the harmonics polarized parallel to the incoming field are 
considered @- 
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Figure 1: (Color online) The spectral amplitudes Q, 2 \a(E p )\ — 
y/S{n)/\d\ (a) and transition dipoles |{0|i|p>| 2 (b) of ID H% 
(the red lines) and its reference atom (the black lines) with 
the continuum |p) obtained by the numerical method. 

According to Ref.[9|], the spectral amplitude a{E p ) is 
\a(E p )\ = ^S(n)/2irn 2 \d\, (1) 

where d(p) = (p|r|0) is the transition dipole matrix ele- 
ment of the atom or molecule between the ground state 
|0) and the continuum state |p), S(£l) is the radiated 
harmonic signal, fl is the photon energy, and E p = p 2 /2 
is the electronic kinetic energy with Q = E p + I p . I p is 
the ionization potential of the atom or molecule. 

Fig. 1(a) plots the spectral amplitudes |a(-E p )| of ID 
with Z=1,R=2 a.u. and its reference atom with 
Z=0.9, R=Q a.u., where |a(_E p )| is calculated from the 
numerical spectrum data by dividing by the transition 
matrix elements |d| = |(0|x|p)| with the |p) obtained 
by numerically diagonalizing the field-free Hamiltonian 
Hq = p 2 /2 + V(x). The corresponding transition dipoles 
used in Fig. 1(a) are shown in Fig. 1(b) and will be 
discussed in detail bellow. The amplitude for in Fig. 
1(a) is consistent with that of its reference atom in the 
high energy region but different in the low energy region, 
because it shows a pronounced peak structure, i.e., the 
amplitude of the peaks is one-order of magnitude higher. 

We now extend these calculations to 2D. Fig. 2 plots 
the spectral amplitudes of 2D models of H^" for varied 
molecular parameters and field intensities. The spec- 
tral amplitudes |a(_E p )| are calculated from the numer- 
ical spectrum data by dividing by the transition matrix 
elements in the length gauge |di en | = |(0|r|p)| with |0) 
obtained by the LC AO-MO approximation, and |p) ob- 
tained by the plane wav e ap proximation and the disper- 
sion relation fl = p 2 /2^~§. For 9 = 90°, the spectral 
amplitudes of the models of and their corresponding 
reference atoms, which have been shifted vertically to 
compensate for differing the overall efficiency, are analo- 
gous (see Fig. 2(a) and (b)). However, for other angles, 
the most remarkable feature exhibited in Fig. 2(c)-2(f) 
is the peak structure at the electron kinetic energy that 
shifts rightward as the orientation angle increases. Ex- 
cept for these peaks, the spectral amplitudes of the other 
energy regimes are analogous. In addition, for the same 
orientation angle but different internuclear distance, the 
location of the peaks is shifted. For example, at 9 = 50°, 
the peak is at E p = 1.85 a.u. in Fig. 2(c), but it shifts 
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Figure 2: (Color online) The spectral amplitudes S7 2 |a(_E p )| = 
■\/S(Q.)/\di en \ of 2D for varied molecular parameters and 
field intensities (indicated in each subpanel) calculated in the 
length gauge with the dispersion relation fi = p 2 /2. In Fig. 
2(a) and 2(b), we also plot the |a(_E p )| of those corresponding 
reference atoms for comparison. 

to E p ~ 3 a.u. in Fig. 2(d). For the same internuclear 
distance, at a fixed orientation angle the location of the 
peaks is almost unchanged when laser intensity varies, as 
revealed by Fig. 2(c) with / = 8.5 x l0 14 W/cm 2 and 2(e) 
with I = 5 x \Q 1A W/cm 2 . These calculations on ID and 
2D for diverse internuclear distance and laser inten- 
sities show that the spectral amplitudes of the molecules 
are usually dependent on the molecular orientation. 

Next, we explain the mechanism behind these phenom- 
ena according to the Lewenstein modelpjj], where the 
time-dependent wave functions can be expanded under 
the strong-field approximation, \tp(t)) = e lIpt [a(t)\0) + 
J dpcp(i)|p)], where a(t) w 1 is the ground-state ampli- 
tude, and c p (t) are the amplitudes of the corresponding 
continuum states, and can be written in the closed form 

cp(i) = i [ dt'E(t') ■ d(p + A(t) - A(O) 

Jo (2) 

xe -*/ t '[(P+A(t)-A(t")) 2 /2+/ P ]rft" ; 

where A(f) is the vector potential of the laser 
field E(f). The time dependent dipole moment 
can be written as D(t) = (ipo(r,t)\r\ip c (r,t)), where 
ipo{ r , t)=e lIpt \Q) denotes the initial electronic state and 
^ c (r,t)=e iZ f* /dpcp(t)|p) denotes the REWP. Eq. 2 
shows that the time-dependent continuum amplitude 
c p (i) is not only dependent on the particular ionization 
potential I p in the exponent of the integrand, but also 
depends on the ground state wave function through the 
transition dipole d(p) in the prefactor of the integrand. 
Thus, the continuum amplitudes of the molecule and its 
reference atom sharing the same ionization potential I p 
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Figure 3: (Color online) The transition dipoles |{0|r|p}| 2 
along the laser polarization direction of 2D with different 
R and orientation angles 9 and their corresponding reference 
atoms, calculated in the length gauge with the dispersion re- 
lation Q = p 2 /2. 



should have the same time-dependent phase, but their 
transition dipoles d(p) differ. 

The above theoretical analysis has been checked by 
numerical simulation. For the ID case, Fig. 1(b) clearly 
demonstrates that the transition dipole for H is different 
from that for its reference atom. The transition dipole 
for shows a deep hollow (the amplitude of the hol- 
low may be up to four order of magnitude lower) around 
E p = 0.33 a.u., but is analogous with its reference atom 
in a broad higher energy region (E p = 1.5 — > 6 a.u.). 
The ionization potential of Hj is similar to its reference 
atom, and their transition dipoles also are analogous in 
a broad high energy region. Therefor, we expect that 
their spectral amplitudes a(E p ), which are determined 
by c p (t), also should be analogous in the high energy re- 
gion, but different in the low energy region, as shown by 
Fig. 1(a). In addition, the location of the hollow in Fig. 
1(b) corresponds to that of the peak in Fig. 1(a). 

Fig. 3 shows the 2D transition dipoles used to pro- 
duce Fig. 2. The transition dipoles for the reference 
atoms (the black curves in each subpanel of Fig. 3) are 
analogous, but the transition dipoles for the models of 
with different R show strong orientation dependence. 
The deep hollow moves as the angle 9 changes. But, at 
9 = 90°, the transition dipoles for the models of Hj" are 
similar to that of the reference atoms in the energy re- 
gion E p =0 — > 7 a.u.. Thus we expect that, for 9 = 90°, 
the spectral amplitude of H^ should be identical with its 
reference atom, but for 9 ^ 90°, they will show a large 
difference in the energy region where the deep hollow in 
the transition dipole of H^ appears. This theory is veri- 
fied by the plot in Fig. 2. 

The ID and 2D numerical results show that intrinsic 
properties of the molecules already enter the REWP. The 
deep hollows in the transition dipoles of H^~ with differ- 
ent 9 can be read from the interference effect as predicted 
by Muth-Bohm et al^jfy. Under the LCAO-MO approx- 
imation, the ground state of the diatomic molecule is 
10} = E l «#o(r-R/2))+6 J |0 o (r-R/2)), where |<fo(r)> 
is the atomic ground state and and bi are the normal- 



ization factors with at = bi for symmetric mixing and 
a,; = —bi for antisymmetric mixing. For symmetric mix- 
ing, e.g., Hj, the transition dipole can be written as 

d(p) = 2cos(p • R/2)^a;s(p) -HRsin(p • R/2) ^ 0^.(3) 

i i 

Heres(p)- /[e^r^Cr)]*, <£o(p) = /[e"^'^^)]*. 
The cos(p • R/2) of the first term in Eq. 3 represents the 
interference between the two cores of the molecules [THjj . 
and the second term is proportional to the internuclear 
distance R in Eq. 3, leading to the breakdown of trans- 
lation invariance. In our calculations, the second term 
is omitted according to Ref . [It], [H, Ell • Thus, the tran- 
sition dipole of the molecule d(p) is different from that 
of the reference atom s(p) due to the interference effect 
from the term cos(p • R/2). The interference that mani- 
fests itself as the hollows in the spectra of the transition 
dipole or the peaks in the HHG spectral amplitudes for 
the molecules, depends on the alignment angle 9 and can 
be located by pRcos(9) — (2n + l)7r where n is an inte- 
ger. For R = 2 a.u. and 9 = 0°, the formula with the 
dispersion relation fl — p 2 /2 predicts the hollow location 
of E p = 0.12 a.u., which is close to our numerical result 
of E p = 0.35 a.u. in Fig. 1(b). If the dispersion relation 
ft = p 2 /2 + is adopted, the location is estimated at 
E p = 1.23 a.u., far away from the numerical result. We 
have calculated the dependence of the minima positions 
(in energy E p ) on molecular internuclear distance R in 
the ID case. It shows that for the larger R (i.e., R > 3 
a.u). and larger momenta p (i.e., p > 3 a.u.), where the 
LCAO-MO approximation for the ground state |0) and 
the plane wave approximation for the continuum |p) are 
expected to be more applicable, the above simple formula 
gives a good estimation on the location of the minima, 
i.e., the deviation is less than 1/10. 

In addition, the transition dipoles for the reference 
atoms with different ionization potentials, as shown in 
Fig. 3(a) and 3(b), are almost identical. We attribute 
the difference in their HHG spectral amplitudes revealed 
by comparing Fig. 2(a) and 2(b) to the difference 
in their ionization potentials. The dependence of the 
spectral amplitude a(E p ) on the ionization potential is 
also evident in the current theories of tunnel ionization 
(ADK[2(]], Yudin and Ivanov 21 . MO-ADKQ), where 
the ionization potential explicitly appears in the expo- 
nent of the formula determining the cycle dependence of 
the ionization rate. Thus, we expect that the time de- 
pendence of the continuum wave packet also relies on the 
ionization potential. Moreover, the present experiment 
for the orientation dependence of ionization rate of CO2 
shows a pronounced deviation from the prediction by the 
MO-ADK theory. The interference caused by the sepa- 
rated cores is not well described by the simple theory (23j. 
Thus, we expect that the dependence of the continuum 
wave packet on the HOMO is also not well described 
by the model. Based on this analysis, we conclude that 
due to the interference effect, the ionization and acceler- 
ation processes in the three-step model[24| also play an 



4 



°-5 



(a) N 2 


— • — 6=10 


— • — 0-1 ° C 


— ► — 0=30° 


(b) N 2 -^6=30' 




— « — 6=40° 


— * — G=40 c 




— t— 0=50° 


— •— 9=50 c 




^r— 0=80° 




y tj 










^ \7 




— •—8=10° 


— e=10 c 


(c) co 2 


— >— 0=30° 


(d) CO z -«~o=30 [ 




— <— 6=40° 


- e=40 c 




— 6=50° 
■ ' "=«"" 


— ♦— e=50' 
* 0=80' 


ft? YvvS 

A ] ♦ \A v 

T l 













■2 ^ 

"a 

2 f 
1 


-1 



3 6 9 1 2 

Electron Kinetic Energy E {a.u.) 



Figure 4: (Color online) The transition dipoles |{0|r|p}| 2 
along the laser polarization direction ((a) and (c)) and spec- 
tral amplitudes Q 2 \a(E p )\ = ^/S(n)/\d len \ ((b) and (d)) of 
N2 and CO2 with different orientation angles 9 calculated 



p /2. The laser intensity is 
1A W/cm 2 and wave length A = 800 nm. The 
(2pa g ) and CO2 (2pir g ) are expressed in the 
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important role in HHG from molecules. In fact, in our 
calculations, the hollows in Fig. 3 for molecules are re- 
sponsible for the interference minima in harmonic spectra 
predicted by Lein et al\j^. 

This analysis demonstrates that the spectral ampli- 
tude of the molecular REWP is different from its ref- 
erence atom due to the interference. However, why 
can the molecular orbital tomography experiment recon- 
struct the HOMO of N2 by its reference atom Ar [9|? In 
Fig. 4, we plot the transition dipoles |d; ell | 2 = |(0|r|p)| 2 
along the laser polarization direction ((a) and (c)) and 
the spectral amplitudes |a(-E p )| ((b) and (d)) of N 2 and 
CO2 with different orientation angles 8. The correspond- 
ing harmonic spectra are obtained by the Lewenstein 
model with the dispersion relation SI = p 2 /20, [l4|. For 



N2 with R = 2.079 a.u., the interference hollows at differ- 
ent angles 9 should all appear in the high energy region, 
i.e., far away from the plateau regime. The curves of the 
transition dipoles in Fig. 4(a) are analogous to each other 
in the low energy region (E p < 2 a.u.). Accordingly, the 
spectral amplitudes in Fig. 4(b) also are analogous to 
each other in the same energy region within a vertical 
scaling factor. However, for CO2 with R = 4.38 a.u., the 
curves of the transition dipoles in Fig. 4(c) show a large 
difference in the low energy region. The curves of the 
spectral amplitudes in Fig. 4(d) also are obviously dif- 
ferent, showing a sharp peak shifting toward the higher 
energy region as the angle 8 increases. We thus expect 
that the spectral amplitude for N2 in the low energy re- 
gion is largely independent of molecular orientation [?|. 
This conclusion is not applicable to CO2, for which the 
interference hollows for small orientation angles 8 should 
appear in the lower energy region corresponding to lower 
harmonic orders in the plateau region 0, Q. The spec- 
tral amplitude of CO2 in that energy region is largely 
related to its orientation. Therefore, for some molecules, 
e.g., CO2, the interference effect should be included in 
the orbital tomography experiment to obtain an accu- 
rate HOMO wave function. 

In summary, from our ID and 2D numerical calculation 
and analytical deduction, we find that the assumption of 
the tomography method, i.e., that the recolliding electron 
wave packet of a reference atom may be used to extract 
the HOMO from high order harmonic spectra, may not 
be generalized to other molecules such as CO2. For these 
molecules, the molecular properties already enter the rec- 
olliding electron wave packet so the spectral amplitude of 
the molecular REWP shows a large difference from the 
reference atom and a strong orientation dependence in 
some energy regions due to the interference between the 
atomic cores of the molecule. 
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